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483. 


ON A PAIR OF DIFFERENTIAL EQUATIONS IN THE LUNAR 
THEORY. 


[From the Monthly Notices of the Royal Astronomical Society, vol. xxxu. (1871—12), 
pp. 201—206.] 


I CONSIDER the differential equations 


4 <A ; (a) +5 +—=km’p {44+ 3 cos (2v — 2mt)}, 
= G T) =jmp? { — è sin (2v — 2mt)}, 


which when j=k=1 give the following equations in the lunar theory (D=t— mt): 
1 a 
5 Lt mt — BEB mt yg me — FB me — 4089 mt — BARRY me — IRE m 


+ cos 2D [m + 42 m + 131 ms + 383 m + 510.5 PEP nah 


080 
T mé + 2737 54 8 6 55483 7 4 2389416723 m8 3523 283 m9 
+ cos 4D [$ mt + m + 18552 ms + e mw + 238 sE4agTd2 m hide m°] 
219 6 9 ¢ no 
+ cos 6D[ 355 m+ 151339 m7 + 29887443 m8 4 98978623957 19] 
+ cos 8D [2704 m8 + 70 33633 m9), 
or as far as m’, 


hee > + $35 mt + 83 ms 2775 mê 87 m7 
p=1 Em + F385 Mm +S ny + APP mo + 4787 m 


— E TG PE TESIR 45 09 m5 — 485178 m6 — 24487949 47 
+ cos 2D[— m- 42 m is m TOP m — fh dre mS — 24385852 n] 
— 3 4 y Pre 
+ cos 4D [— 3 mt — 4211 m” — 74069 më — 1749179 m7] 
_ go 
+ cos 6D [— $35 me — 1261.93 m7], 
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eee ae 
(; is given by M. Delaunay only as far as m’, the additional terms of „i and expression 


for p were kindly communicated to me by Prof. Adams) ; and 
v=t 
+ sin 2D (41 m? + 52 m + 523 mt + 2855 m4 8304449 9778 
8051418161 m?) 
+ sin 4D (394 mt + $49 m” + SAIR23 ms 4 8136: 
+ sin 6D (3748 m° + 64511 m7) 
(Delaunay, t. 11. pp. 815, 836, 845). 
To integrate the original equations write 
P=l +p +p + 
v=t kU tt Fo, 


where the suffixes indicate the degrees in the coefficients k, j conjointly: the equations 


for Pn, Vn take the form 
d dpr ae 
dia. a 


d (Fe 


2% 
=e s+ CAEN 


dt \ dt 


where Va, Un, Pn, Qn do not contain p, or v,. From the second equation we have 


dv, 
dt 


+ 2p, + D = Pu 


+ 2Pn + U,=9,+ | P, dt, 


where ©, is a constant of integration, the integral [Pn dt containing only periodic 


terms; and then adding twice this to the first equation we have 


d dpn 


ae at + Pn + Vu +2U,=2% +Q,+2/ Prat 


; Dias oe i ; dv 
which determines p,; and substituting its value in the other equation we have a} 


i : dv ; 
and thence v,; the constant Q, is determined so that —” may contain no constant 


dt 


term. We have 


Veh | | U,= m 
— 2p, A + pip: — 4p, | 
&e. | &e. 
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Qı = km? ($ + $ cos 2D), P, = jm? (— 3 sin 2D), 
Q: = km? {3v sin 2D + p, (4 + $ cos 2D)}, P, = jm? (— 3v, cos D — 3p, sin 2D), 
Q; = km? {— 3v, sin 2D — 3v,2 cos 2D P, = jm? {— 3v, cos 2D + 3v? sin 2D 
+ pv, .3 sin 2D — 6p,v, cos 2D 
+ ps (4+ $ cos 2D)}, + (2p.+ p°). — $sin 2D}, 
&e. &e. 


In particular attending to the values of P,, Q, the equations for pı, v are in their 
original form 


d d dv 

ai 4 3p: +2 ea = km? (4 + $ cos 2D), 
d de, ei ee 

dt (a = 2p.) = jm’ ( — į sin 2D), 


whence in the transformed form they are 


dv, 


di + 2p, = 0, + “n° 2D, 


4 i j= 
and 
ga + pı = 20, + km?($ + § cos 2D) +5 bjm I me” 2D. 
, it a W0 
Thus the constant term of pı is 20, + km’, giving in J ê constant term — 80, = km? 
this must vanish, or we have 0,=—4km?; and the equations thus become 
dv 


3jm? 
ae o a kat TET 8 2D, 


4 


lt + h= tkm ++ (3 km? + 5 bgm m) cos 2D, 


and then completing the integration 


= 3 km? ao 3 m2 
os sb 2 2d 2D 
pe temel E ta -a T] cos 2D, 
Si 3 km? gym? (T= ere 2D, 
ai: lama 8m + 4m) * (1 — mY (3 -8m + 4m’) aM 


which are the accurate values of p, and 2. 


Expanding as far as mê we have 
=k(=1m*)+cos2D{ k(=4 m= 4 m= 2% m YP m YR m) 


s 2 5 4. 575 5 xu 6 
+j(— mim 88 mt a m— SEB m), 


r é P 1; i: = boot 6 
which for j= k is =k(— m?—412 m= pim 828 mi r m’), 
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and 
Y= sn2D{ k( 4 m?+4L m+ 88 m+ 525 m+ 398) m’) 
1 2 18 54 
. t ‘ 6 
Te. i m+ i 2g ym +45 BRO mt + ie m? + 22231 m')}, 
which for j= k is = k 4. m+s 5y Lm + 222 mi S96 m + Ag mi). 


I have, not in general, but for the value j= k, calculated p, and v, as far as am: 
I have not made the calculation for p, and 2, but their values may be deduced from 
the foregoing values of p, v; the final expressions (when j= k) of p, =1 + pit P+ pP: + 


and v, =t +u +U + U... are 


p=1 

+k (—im ) 
T Jajmi+ $p m+ A m) 
+1? ( | — 4984 mì) 

+cos2D{ k (— m?—12 m — 1382 mt — 828 m — 552% m?) 
+ h?( 2 m+ 31 mt 322 m’) 
+t — 3884 m) 

+cos4D{ k( — } m — 1P m” — Bape? m.) 
+k ( + o m’) 

+cos6D{ k ( — #2 m’), 

and 
v=t 

+sin2D { k (Adm? +m 222 mt+ 398 mi + 44975 m’) 
+ k( — ¢m—= 44 m- #4 w) 
+ — $f m) 

+sin4D{ k?( 201 ms + S48 m + 585963 m8) 
+ ke ( — $3 m’)} 

+sn6D{ k( + 2215 m’); 


which for k=1 agree with the foregoing formule (verifying them as far as më); the 
present formule exhibit the manner in which the expressions depend on the several 
powers of the disturbing force. 
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